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ABSTRACT 

< 

O Context. The Milky Way (MW) bulge is a fundamental Galactic component for understanding the formation and evolution of galaxies, 
in particular our own. The ESO Public Survey VISTA Variables in the Via Lactea is a deep near-IR survey mapping the Galactic bulge 
and southern plane. Particularly for the bulge area, VVV is covering ~ 315 deg 2 . Data taken during 2010 and 201 1 covered the entire 
Q j bulge area in the JHK S bands. 

Aims. We used VVV data for the whole bulge area as a single and homogeneous data set to build for the first time a single colour- 
O ' magnitude diagram (CMD) for the entire Galactic bulge. 

Methods. Photometric data in the JHK S bands were combined to produce a single and huge data set containing 173, 150, 467 sources 
in the three bands, for the — 315 deg 2 covered by VVV in the bulge. Selecting only the data points flagged as stellar, the total number 
of sources is 84, 095, 284. 

Results. We built the largest colour-magnitude diagrams published up to date, containing 173.1+ million sources for all data points, 
and more than 84.0 million sources accounting for the stellar sources only. The CMD has a complex shape, mostly owing to the 
complexity of the stellar population and the effects of extinction and reddening towards the Galactic centre. The red clump (RC) 

giants are seen double in magnitude at b 8° - 10°, while in the inner part (b 3°) they appear to be spreading in colour, or even 

splitting into a secondary peak. Stellar population models show the predominance of main-sequence and giant stars. The analysis of 
the outermost bulge area reveals a well-defined sequence of late K and M dwarfs, seen at (J — K s ) ~ 0.7 - 0.9 mag and K s > 14 mag. 
Conclusions. The interpretation of the CMD yields important information about the MW bulge, showing the fingerprint of its structure 
and content. We report a well-defined red dwarf sequence in the outermost bulge, which is important for the planetary transit searches 
of VVV. The double RC in magnitude seen in the outer bulge is the signature of the X-shaped MW bulge, while the spreading of the 
RC in colour, and even its splitting into a secondary peak, are caused by reddening effects. The region around the Galactic centre is 
£N) ■ harder to interpret because it is strongly affected by reddening and extinction. 
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5—1 " Introduction tre are difficult because they are affected b y crowding and ex- 

! ^ , , „ , ~ mTX . r , , ^ , tinction. Optical su rveys such a s MACHO (lAlcock et al.l ll999) 

The bulge of the Milky Way (MW) is a fundamental Galactic and QGLE (Udalsk i etal.ll99l are highly affected by extinction 

component for understanding the formation and evolution not while past near . infrared (near -IR) surveys such as 2MASS are 

only of our Galaxy, but of galaxies in general. In the MW limited tQ brfght sourcefj {K „ U3 S krutskie et al. 2006), 

the stars can be individually resolved down to faint magm- and do nQt allow a complete view of the bulge popu i ations . 
tudes, allowing a complete understanding of its structure, age, 

kinematics, and chemical composition. However, because the The new ESO Public Surve y VISTA Variables in the 
bulge and plane concentrate not only most of the stars in the Vl ' a L^tea (VVV) is a deep near-IR survey mapping 562 
MW, but also gas and dust, observations of the Galactic cen- square degrees in th e southern plane and bulge of our Galaxy 
(IMinniti et al.ll2010 ). Particularly for the bulge area, VVV is 

Send offprint requests to: R. K. Saito: rsaito@astro.puc.cl covering -315 deg 2 . The VVV survey observes in five near- 

* Based on observations taken within the ESO VISTA Public Survey *R passbands (ZYJHK S ), where extinction effects are lower than 

VVV, Programme ID 179.B-2002. in the optical wavelengths, and it is much deeper than previ- 

** The VW survey data are available thr ough the ESO archive ous near-IR surveys. The 5a limiting magnitude of the VVV 

http://www.eso.org/sci/archive.html data using aperture photometry is K s ~ 18.5 mag in clean 
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Fig. 1. Left-hand panels: total number of sources present in the J, H, and K s singleband and in the JHK S multi-band catalogue 
for each VVV bulge tile. The tiles are numbered from b201 to b396, starting at the bottom-left corner of the VVV bulge area 
(/, b — -10, -10). The Galactic latitude for the tile centres is shown in the figure, while dashed lines mark the edge of each line of 
tiles across the bulge. The right-hand panel presents a histogram with the ratio between the number of sources in the multi-band 
catalogue in relation to its least numerous single-band catalogue. 



fields, which allows one to monitor bright sources such as RR 
Lyrae and clump giants stars along the whole Galactic bulge 
and plane, solar-type stars at the Galactic centre distances, and 
even faint dwarf stars a few kpc away. The VVV data reach the 
main-sequence turn-off (MSTO) even in intermediat e- to high- 
extinction areas (E(B - V) < 4.0. lMinniti et al.ll2010h . 

The VVV photometry for the whole bulge area was com- 
bined into a single, huge, and homogeneous data set, allowing 
us to build for the first time a single colour-magnitude diagram 
for the entire Galactic bulge (~ 315 deg 2 ). Here we present the 
84-million star VVV colour-magnitude diagram (CMD) for the 
Galactic bulge. This is the largest CMD ever published for a 
large homogeneous data set, a significant step forward since the 
publicatio n of the 9-million star CMD of the Large Magellanic 
Cloud dAlcock et al.l |2000T) . We discuss the differences in the 
morphology caused by crowding, extinction and sky brightness. 
A forthcoming paper will present a similar analysis for the pop- 
ulations across the VVV Galactic plane area. 

2. Observational data and catalogues 

The VVV is an ESO Public Survey scannin g the MW bulge an d 
plane with the 4-m class VISTA Telescope (IMinnifi et al .11201 Oh . 
The total observed area is about 562 deg 2 , within -10.0° < 
I < +10.5° and -10.3° < b < +5.1° in the bulge, and within 
294.7° < I < 350.0° and -2.25° < b < +2.25° in the plane. The 
VVV survey observes in five passbands, namely Z (0.87 /im), Y 
(1.02 yum), J (1.25 fim), H (1.64 /mi), and K s (2.14 jjm) bands, 
and also conducts a variability campaign in /T s -band only, with 
~100 pointings spanning five years (2010-2014). 

Each unit of VISTA observations is called a (filled) "tile", 
consisting of six individual (unfilled) pointings (or "pawprints") 
and covers a 1 .64 deg 2 field of view. To fill up the VVV area, a 
total of 348 tiles are used, with 196 tiles covering the bulge (a 
14 x 14 grid) and 152 for the Galactic plane (a 4 x 38 grid). 

The schedule for the first year (2010) comprised the ob- 
servation of the whole VVV plane and bulge areas in the five 
passbands, as well as five complementary epochs in /T,-band. 
However, delays during the campaign prevented the completion 
of the planned observations, with some of the first year data be- 
ing observed during 201 1. 



The delays mostly affected the Z and Y observations and the 
variability campaign, while the JHK S observations were fully 
completed during the first semester of 201 1 . From the 196 point- 
ings needed to cover the total bulge area, 180 were completed 
and made publicly available with the VVV Data Release 1 (DR1) 
on August 2011 (a compre hensive description of this dataset is 
given by ISaito et al.ll2012h . The present work is based on J, H 
and K s photometry, from the VVV DR1 and DR2 data (in prepa- 
ration), with complete bulge coverage in the three bands. 

Photometric catalogues for the VVV images are provided by 
the Cambridge Astronomical Survey Unit (CASUfl The cata- 
logues contain the positions, fluxes, and some shape measure- 
ments obtained from different apertures, with a flag indicating 
the most probable morphological classification. In particular, we 
note that "— 1" is used to denote the best-quality photometry of 
stellar objects. Some other flags are "-2" (borderline stellar), 
"0" (noise), (non-stellar objects), "-7" (sources containing 
bad pixels), and "-9" (saturated sources). Our entire analysis 
was performed using the apermag3 aperture fluxes, which are 
used in the CASU catalogues and in the VISTA Science Archive 
(VSAfl as the default values to represent the flux. However, we 
note that the CASU pipeline starts to measure the positions and 
fluxes with the aperraagl (1" diameter) and then successively 
increases the aperture by a factor of V2 in diameter. For highly 
crowded fields as in the inner Galactic bulge, using apermagl 
and apermag2 can be more suitable than apermag3. 

The VVV data are in the natural VISTA Vegamag system, 
with the photometric calibration in JHK S performed using the 
VISTA magnitudes of unsaturated 2MASS stars present in the 
images. A final tile catalogue comprises the photometry for 
an area of ~1.64 deg 2 in the sky, with the total number of 
point sources found in each single-band catalogue ranging from 
408, 858 at the outermost bulge region (tile b202, Zf s -band) to 
1, 597, 046 sources per tile at the innermost Galactic centre (tile 
b246, /-band). The /T s -band catalogues are shallower than J 
and H at high Galactic latitudes, while near the Galactic plane 
(\b\ < 3 deg) the three catalogues present a similar number of 
sources for most tiles (see Fig. [TJ. On the Galactic plane the 



1 http : //casu . ast . cam . ac . uk/vi stasp/ 

2 http : //horus . roe . ac . uk/ vsa/ index . html 
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Fig. 2. Density plot in logarithmic scale showing the VVV bulge area. The map was made using the multi-band JHK S CASU vl.l 
and vl.2 catalogues for point sources brighter than K s = 16.5 mag (see text). Crowded areas appear in yellow, while less populated 
regions as well as highly extincted areas are shown in blue. The overlapping regions between the tiles are highlighted since the point 
sources are accounted for twice. Dashed lines mark the 15 areas selected to build the colour-magnitude diagrams shown in Fig. [6] 
(see Section 4). The source density, in units of 10 s sources deg -2 , is indicated in the horizontal bar at the top. 



high extinction mostly affects the J- and //-band, which ren- 
ders the ZT s -band catalogue the one with the most entries among 
the three bands. 

The comparison between the source counts of different tile 
catalogues depends strongly, but not exclusively, on their posi- 
tion with respect to the stellar crowding, i.e., the detection rate 
for tiles at high Galactic latitude will be significantly higher than 
for a tile in the Galactic plane. This effect becomes even more 
important for the aperture-photometry-based source catalogues 
used in this study. Artificial star experiments carried out on sub- 
sections of two different tiles in t he bulge area (b204 and b314, 
respectively, see lSaito et al]l2012l) show that the single-band tile 
catalogues are 80% complete down to K s ~ 17.0 in the outer- 
most bulge region (b204), but in the innermost bulge (b314) the 
same percentage is reached at K s ~ 15.5. 

The single-band CASU catalogues for each tile were 
matched by VVV team members, with most of sources matched 
within < 75, as expected fr om the astrometric accuracy of the 
VVV data dSaito et alj|20 12). For a couple of tiles the elliptic- 
ity sometimes varies between the filters, hence sources appear 
as high-quality stars in one filter but are rejected in another, or 
appear with a distinct morphological classification. Therefore, 
some sources are lost during the matching procedure. For most 
tiles the multi-band catalogues contain more than 85% of of the 
total sources found in the least numerous singleband catalogue 
used for the matching (see Fig.Q]). The performance is maximum 



at high Galactic latitude, with the multi-band catalogues reach- 
ing up to ~96% of successful matches. On the other hand, the 
increasing source confusion and extinction towards the Galactic 
centre reduces the number of matches in the innermost region 
(\b\ < 1 deg), with a few tiles showing ~75-80% of the to- 
tal sources present in the least numerous singleband catalogue. 
Some larger "steps" seen between adjacent tiles in Fig. Q] are 
caused not only by differences in the quality and deepness of 
data, but also by the obscured regions at lower Galactic latitudes. 
There are 173, 150,467 sources in all multi-band catalogues for 
the 196 bulge tiles. 

Figure [2] shows the VVV bulge area as a density map corre- 
sponding to the mosaic of all multi-band catalogues. To equalize 
the increasing sky brightness towards the Galactic centre that 
is generated by the contribution of the underlying, unresolved 
faint stars, a cut in K s = 16.5 mag was applied in the cata- 
logues, with a total of 128,660,076 sources remaining in the 
map. High-density areas appear in yellow, while less populated 
regions and high-extinction areas are shown in blue. The tiling 
pattern produces overlapping regions between the tiles. These 
are highlighted in the map, since the point sources in these re- 
gions are counted twice. The overlapping area has been used by 
the VVV team to check the photometric and astrometric accu- 
racy, but it was not been coadded. 

In a few cases the colour pattern (or source density) changes 
between adjacent tiles, because some sources in lower quality 
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Fig. 3. K s x{J - colour-magnitude diagram for the VVV bulge area, as shown in Fig. [2] In the left-hand panel is the CMD for all 
point sources found in the matched catalogues, in a total of 173, 150,467 sources. The right-hand panel shows only point sources 
flagged as "stellar" in the three J,H, and K s catalogues, according with the CASU photometry. The total number of stellar sources is 
84, 095, 284. Contour lines mark density levels in steps of 5% from the maximum density. The red dening vector associa ted with an 
extinction of E(B —V) — 2, based on the relative extinctions of the VISTA filters, and assuming the|CardemetaD([l989) extinction 
law, is also shown in both panels. 



Table 1. Effective wavelengths for the VISTA filter set used 
in the VVV obser vations and the relati ve extinction for each 
filter based on the ICardelli et al.1 (1 19891) extinction law (from 
ICatelan et al.ll2011h . 



Filter 




A x /A v 


A X /E(B - V) 


J 


1.254 


0.280 


0.866 


H 


1.646 


0.184 


0.567 


A' s 


2.149 


0.118 


0.364 



images were rejected, as cited above. This effect is mostly seen 
near the Galactic centre, where the detection rate is more affected 
by crowding. Several tiny dark spots seen in the map mark the 
regions with no detections around bright, saturated stars. The 
globular cluster M22 is also detected at /, b — 1-9.9, -7.5. It ap- 
pears as a region with increasing density, which suddenly drops 
to a spot with no detections, due to confusion of the aperture 
photometry in detecting sources in highly populated fields such 
as globular cluster cores. 

3. The 84-million star CMD 

We used the JHK S VVV multiband catalogues of the complete 
bulge as a single data set to build the largest colour-magnitude 
diagrams published to date. The left-hand panel of Fig. |3]shows 
the CMD containing 173, 150,467 sources for the ~ 315 deg 2 
area covered by VVV in the bulge. It contains all sources found 
in the catalogues and it is shown as a density plot with con- 
tour curves marking the same star density levels in steps of 
5% from the maximum intensity. The CMD shows the same 



data as in Fig. [2] but also includes the point sources fainter 
than K s = 16.5 mag, in contrast with the additional selection 
criteria applied in Fig. [2] The brightest stellar sources appear 
at K s ~ 9.9 mag, while the 5<x limiting magnitude depends 
strongly on crowding. In the outer bulge the limiting magnitude 
is K s ~ 18.3 mag, but at the Galactic centre the presence of un- 
derlying, unresolved faint stars contributes to increase the sky 
brightness, which limits the photometry to K s ~ 16.3 mag. A 
small fraction of the sources, flagged as "saturated stars" and 
"noise", appears to be populating the upper and lower limits of 
the CMD, respectively. 

The right-hand panel of Fig. [3] shows the CMD built only 
with the data points flagged as "— 1", stellar sources, in the J, 
H, and K s catalogues, in a total of 84,095,284. Since the stel- 
lar flag comprises only best-photometry data points, structures 
are clearer seen in the right-hand panel of Fig. [3] Moreover, the 
limiting magnitude is shallower, since the faint sources flagged 
as "noise" are not present. Fig. [3] also shows in both panels the 
reddening vector associated with an extinction of E(B — V) =2, 
based on the relative extinctions of the VISTA filters using the 
ICardelli et alJ dl989l) extinction law (see Table 1). However, we 
notice that the reddening law varies significantly in the inner- 
most regions of the MW (e.g., Nishivama et al. 2009). 

Because the CMDs were created by combining all individ- 
ual 196 tile catalogues, they include the overlapping regions be- 
tween the tiles, which comprise ~ 7% of the total point sources. 
Thus, ~ 12.6M of the 173. 1M sources are counted twice, which 
reduces the number of unique sources to ~ 160.4M. Similarly, 
~ 78. 0M of the 84. 1M stellar sources are unique. The overlap- 
ping areas are equally spaced along the bulge in both I and b 
directions, and therefore do not contribute any bias or special 
trend in magnitude or colour. 
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Fig. 4. Top-left panel: K s x (J - K s ) CMD for an area of 1 sq. deg. centred at (/, b) = (0.0, -9.5). This is the outermost region 
in the VVV bulge area, on the Galactic minor axis. The different stellar populations are labelled. Red circles show the colours of 
M-dwarfs for a distance modulus of (m^ s - Mk s ) = 7 mag (see text). The top-central and top-right panels show stellar population 
synthesis models of the Galaxy for the same region, provided by the Besan^on (centre) and Trilegal models (right-hand), including 
photometric errors and an extinction as described in the text (see Section 4). In the Besancon model stars are shown according to 
luminosity class. Main-sequence stars are depicted in orange, subgiants in blue and giants are red dots. Supergiants, bright giants, 
and white dwarfs are also present, but there are only a few (show in light blue, magenta and purple, respectively). The relative 
numbers are given in Table 3. On the other hand, the Trilegal model highlights the Galactic components, with 69.6% of the stars 
belonging the bulge (in red), 18.3% are located in the halo (in blue) and 10.6% are thin-disk stars (in orange). A small fraction of 
1.5% are thick-disk stars (green dots). The bottom-left panel shows the luminosity function for data and model. The bottom-right 
panel shows the luminosity function of each luminosity class in the Besan§on data. The colour code is the same as used in the 
CMDs. 



The following description, and the next sections discuss the 
high-quality CMD that presents the stellar sources. 

The VVV Bulge CMD has a complex shape, mostly be- 
cause of the complexity of the stellar population, seen at dif- 
ferent ranges of magnitudes and colours, and the effects of ex- 
tinction and reddening towards the Galactic centre. Reddening 
affects the bulge populations at different levels, spanning from 
E(B - V) < 0.2 atfr ~ -10°, and up t o E(B - V) ~ 10 at the the 
Galactic centre dGonzalez et alj|2012h . Moreover, the reddening 
law changes in the innermost bulge regions, as cited above. 

In particular, the red giant branch (RGB) is seen to be 
broader, with the main peak of the red clump (RC) giants at 
(J - K S ),K S ~ 1.05, 13.20. The RC giants have structures rang- 
ing in colour and magnitude, with a secondary peak in colour 
at (J - K s ), K s ~ 0.85, 13.20, and an elongated structure in 
magnitude with AK S ~ 1 mag. Those are caused mostly by the 
X-shaped bulge of the MW, producing the structure in magni- 
tude, and by differential extinction, resulting in multiple struc- 
tures in colour (see Section 5.1). The colour of the RGB stars is 
also affected by changes in the m etallicity, which sh ows a gra- 
dient along the bulge minor axis (IZoccalietal.1 12008). However, 
the effect is minimized in the RC, since the c olour of RC giants 
is les s affected by age and metallicity (e.g., ISalaris & Girardil 
2002). 



4. Comparison with models 

In this section we compare our data with stellar popu lation syn- 
thesis models provided by the Be sancorPI (R obin et al.l2003l) and 
the Trilegafl dGirardi et al.ll2005l) codes. 

For the comparison we selected three regions of 1 x 1 along 
the bulge. The first region is at high Galactic latitude, centred at 
(/, b) = (0.0, -9.5). This is the outermost region of the VVV 
bulge area, and is the least crowded as well as least affected by 
reddening on the bulge minor axis, with A^ s = 0.023 ± 0.010 
(see Table 2). Owing to these conditions its stellar population is 
well defined, allowing a better comparison with models. 

The other two selected regions are centred at (I, b) = 
(0.0, -6.0) and (/, b) = (0.0, -3.0). These regions are highly af- 
fected by reddening and crowding, and are useful for testing the 
differences in the morphology of the CMDs. The mean values of 
reddening and extinction for each region are given in Table 2. 

To properly account for the differential reddening across the 
1 sq. deg. regions, we divided each region into 16 small areas 
of 15' x 15' side. The mean Ak s , E(J - K s ), and respective 
Ay were calculated for each area using the reddening maps of 



3 http://model.obs-besancon.fr/ 

4 http : //stev . oapd . inaf . it/cgi-bin/trilegal_ 1 . 5 
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Fig. 5. Comparison between VVV and synthetic data provided by the Besancon model for the areas at (I, b) — (0.0, -6.0) (top), and 
(/, b) — (0.0, -3.0) (bottom). The notation is similar to that presented in Fig. [4] 



Table 2. Mean values for the extinction and reddenin g for each 
1 sq. deg. region pr ovided by iGonzalez et alJ (1201 21) using the 
ICardeili et aild 19891) extinction law. 



(l,b) [deg] 


A K s 


E( J - 


K s ) 


A 


V 


(0.0,-9.5) 
(0.0, -6.0) 
(0.0, -3.0) 


0.023 ± 0.010 
0.110 ±0.014 
0.320 ± 0.052 


0.033 ± 
0.160 ± 
0.464 ± 


0.014 
0.020 
0.075 


0.194 + 
0.932 ± 
2.710 ± 


0.083 
0.119 
0.438 



IGonzalez etal] d2012l) . and the ICardelli et ail d!989l) extinction 
law. These values were used as input while querying the models 
for the synthetic data of each 15' x 15' area. Finally, the data 
generated by the models for each 16 small areas were combined 
into a single data set corresponding to the entire 1 sq. deg. This 
procedure ensures that the spatial distribution of sources across 
the region and the changes in the reddening are well-accounted 
for the synthetic data. 

The Trilegal model is available in the VISTA JHK S colours, 
while the Besancon model uses the JHK system. Although our 
data use the VISTA system, the systems are sufficiently close 
that a comparison is useful for the conclusions. When synthetic 
Besancon models become available in the VISTA system, these 
should be used instead. The synthetic data were collected us- 
ing the default parameters presented in both models, which in- 
clude all luminosity classes belonging the four Galactic popu- 
lations (thin disk, thick disk, halo and bulge) at the given line 
of sight, photometric errors, and extinction as described above. 
While the Besancon data output labels the sources according to 
the luminosity class, the Trilegal model complements the anal- 
ysis because it classifies the sources according to the Galactic 
component. 



The top-left panel of Fig. g] shows the K s x (J - K s ) CMD 
for the 1 sq. deg area at (/, b) = (0.0, -9.5), with labels in- 
dicating the distribution of the stellar content. For quantita- 
tive purposes, the data include not only the stellar sources, but 
all unsaturated sources presented in the JHK S multi-band cat- 
alogue. The CMDs built with Besancon and Trilegal data are 
shown in the top-central and top-right panels of Fig. [4] respec- 
tively. Complementary K s luminosity functions are also shown 
in Fig. |4] In the bottom-left panel we present the luminos- 
ity function of all sources in each model, compared with the 
VVV data. The luminosity functions are quite similar, but the 
Besangon data are most numerous than VVV (+37% sources in 
the range 13 < K s < 17), while Trilegal delivers less stars than 
the real data (-49% in the same K s range). This discrepancy in 
total numbers is expected because models based their calibra- 
tion on data with different levels of completeness in the magni- 
tude range of the VVV data (e.g., 2MASS, DENIS, OGLE, etc). 
However, the difference in the number of sources between VVV 
and the Besancon luminosity function entirely agrees with the 
results of our completeness tests performed on the VVV data 
in this region (see Section 2). The bottom-right panel of Fig. [4] 
presents the luminosity function of the Besancon data separately 
for each luminosity class. 

We notice that the double-peaked RC is not modelled by the 
synthetic data and therefore does not allows a direct compari- 
son. A recent paper by the Besancon team discuss the modeling 
of the double clump by introducin g of a flaring bar in the current 
version of the Besancon code, see Rob in et alJ d2012l) . However, 
since they compared the model with 2MASS data, which are not 
complete at the fainter magnitude of the second clump, addi- 
tional modifications in the model based on deep near-IR data are 
necessary (e.g., VVV and UKIDSS-GPS). 

The Besangon model represents the VVV data well, showing 
the disk and bulge sequences at similar colour and magnitudes, 
as well the similar scattering at lower magnitudes due to the large 
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Table 3. Relative number of each luminosity class in the syn- 
thetic data provided by the Besancon model. The range in 
was defined according to the limiting magnitude of the VVV 
data in each region. 



Luminosity classes 
(relative number [%]) 


QM [deg] 


(0.0, -9.5) 


(0.0,-6.0) 


(0.0,-3.0) 


K s range [mag] 


11.5-18.0 


11.5- 18.0 


11.5-17.0 


Main-sequence 


76.9 


67.7 


17.8 


Subgiants 


16.5 


22.7 


35.9 


Giants 


6.5 


9.5 


45.9 


Supergiants 


0.04 


0.02 


0.05 


White dwarfs 


0.05 


0.06 


0.3 


Bright giants 


<0.001 


0.001 


0.02 



photometric errors. Main-sequence stars and subgiants are the 
most numerous sources, accounting for 76.9% and 16.5% of the 
total number, respectively. Giants account for 6.5%, but domi- 
nate the CMD for K s < 14 mag. A small fraction of luminous 
supergiants, bright giants and white dwarfs is also present. Table 
3 summarises the relative number of each luminosity class in the 
Besangon data. Finally, the top-right panel of Fig. |4] shows the 
CMD for the Trilegal data. The sources are labelled according to 
their Galactic components, which are distributed as follows. The 
bulge population is predominant with 69% of the total sources. 
Halo stars from farther out, behind the bulge, account for 19.6%, 
while the thin disk contributes 10%. Thick disk stars make up 
1.4% of the total. 

The comparison of our data with the models for the region 
at (/, b) = (0.0, -9.5) allows us to conclude that ( f) the MSTO 
highlighted in the VVV CMD is composed of bulge, halo and 
thin-disk stars, with a majority of main-sequence stars, but with 
a small fraction of subgiants. (if) The sequence labelled at the 
left side as main-sequence (MS) stars is composed of subgiants 
and main-sequence stars from the thin disk, (fff) The bulge RGB 
at the right side is composed of giants from the bulge, but also in- 
cludes a small fraction of halo stars behind the bulge, (fv) Finally, 
a sequence seen at the right-most side is comprised of late dwarfs 
belonging to the thin-disk. 

Fig.|5]shows the VVV data compared with the synthetic data 
from the Besangon model for the regions at (/, b) - (0.0, -6.0) 
and (I, b) = (0.0, -3.0). The increasing crowding and reddening 
towards the Galactic centre introduces degeneracies both in mag- 
nitude and colours, that is why at b — -3.0 most of the sequences 
are seen overlaid. The relative number of each luminosity class 
also changes, not only due to intrinsic changes in the stellar con- 
tent, but also because the increasing limiting magnitude of the 
VVV data that contribute to reduce the number of faint sources. 
While the number of main-sequence and subgiants decreases to- 
wards the Galactic centre, the bulge giants dominate the CMD 
at the inner bulge, increasing by a factor ~ 8 from the CMD at 
b = -9.5 (6.5% of the total sources) to the CMD at b = -3.0 
(46%). A more detailed comparison with the models is given in 
the next section. 



5. Slicing the VVV bulge CMD 

To study the bulge CMD in detail we divided the VVV bulge 
observations into 15 smaller areas, each comprising between 8 to 
15 tiles (~ 12 - 22 deg 2 , respectively), as marked in Fig. [2] This 



analysis allows us to check by a quick view the main changes 
caused in the bulge stellar content by the differential reddening 
and extinction towards the Galactic centre. 

Fig. [6] shows the 15 CMDs, with the areas numbered from 
#1 to #15, starting at the bottom-left corner. The mean stellar 
density and the central coordinates for each area are shown in the 
panels. Area #1 contains not only MW stars, but also the core of 
the Sagittarius dwarf galaxy, in a small prop ortion of ~ 1 : 1000 
stars, according with Cseres nies et al.l (|2000). 

The main differences between the CMDs occur at the 
Galactic latitude, while differences in Galactic longitude (which 
spans from -10 < / < +10) at the same b produces smaller 
changes. Thus, the central CMD in each line of panels of Fig. [6] 
corresponding to a given area along the bulge's minor axis, is 
representative of its neighbours at higher longitudes. Fig.|7]high- 
lights the region around the red clump giants for these central 
areas on the minor axis, from the outermost bulge (Area #2) 
to the Galactic centre region (Area #11). The figure shows the 
K s x(J -K s ) CMDs and the (/ - H) x (H - K s ) colour-colour 
diagrams (CCDs) for the given areas. 

The RC of Area #2 is found to be double-peaked, with 
the brightest peak at K s = 12.75 mag and the faintest one at 
K s = 13.45 mag. The peaks appear with similar intensity to each 
other. The CCD of area #2 shows the RC as single-peaked, but 
displaced from the centre of a peanut-shaped structure. This can 
be interpreted as the sum of two closer distributions, as seen in 
the respective CMD. 

The second row from the bottom of Fig. |6]shows the CMDs 
for areas #4, #5 and #6. These areas are moderately affected by 
extinction and reddening, and the main structures appear at sim- 
ilar colours and magnitudes as in the outermost region (areas 
#1 - #3). The RC becomes more prominent in relation to the MS 
and it is seen extended in magnitude. The CMD of Fig.|7]reveals 
a double-peaked structure, but now the farther RC (the fainter 
one, K s = 13.40) is stronger than the closest one (the brighter, 
K s = 12.95), compared with area #2, where the peaks show simi- 
lar intensity. The peaks are closer to each other compared to area 
#2, with AK S = 0.45, compared with AK S = 0.70 for the outer- 
most bulge. The CCD of Area #5 shows a similar behaviour as 
area #2, with a peanut-shaped structure with a single maximum. 

Areas #7, #8 and #9 are affected by reddening and extinction 
at different levels, and show the RC spread over a broader colour 
range. Because the star density increases, the unresolved stars 
contribute to increase the background brightness limit, with the 
faintest point sources ~ 0.5 mag brighter around the minor axis 
(area #8). The red giant branch is seen single in area #7, with the 
RC also single, but spreading along the reddening vector. In area 
#8 the RC, and even the RGB seems to be double in colour, but 
in a closer separation. This separation becomes larger in area #9. 
The inspection of Fig.|7]shows that the RC main peak for area #8 
appears at (/ - K s ) ~ 1.05 over an elongated structure reaching 
(J - K s ) - 0.90, but at similar magnitude (K s ) - 13.15 mag). 
Its CCD shows similar behaviour, with a single peak over an 
elongated peanut structure. 

The Galactic centre is sampled in areas #10, #11 and #12. 
These are harder to interpret, since they are strongly affected 
by reddening and extinction. The sky brightness is even higher 
around the minor axis, with the faintest objects seen at K s ~ 
16.5 mag in area #1 1 . This is the reason for the magnitude colour 
cut applied in Fig. [2] All populations are seen to be much redder, 
with the RC spreading along the direction of the reddening vec- 
tor. The RC has a main peak at (J - K s ) ~ 1 .50, with a difference 
of A(J - ~ 0.8 mag in relation to the outermost region. 
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Fig. 6. CMDs for 15 different areas along the bulge as shown in Fig. [2] The areas were labelled from #01 (bottom-left corner) to 
#15 (top-right). The mean density of stellar sources (sources deg~ 2 ) and the central coordinates (I, b) for each area are shown in the 
panels. 



Areas #13, #14 and #15 are centred at b — +3.54 deg, and 
show a similar behaviour as seen for areas #07-09 (centred at 
b - -3.01 deg) but they are slightly less affected by extinction. 
At the minor axis (area #14) the RC stars are more prominent 
compared to the MSTO, while in the higher latitudes (areas #13 
and #15) both RC and MSTO have similar strength. 



5.1. The red clump giants 

Owing to th e well-defined absolute magnitude and colour (e.g., 
lAlvesll20"ooT) . and the well-understood de pendence on age and 
metallicity (e.g., [Salaris & Girardil l2002h . clump giants have 
been used with su ccess as distance indicators and calibrators in 
extinction studies (IPaczynski & Staneklll998tlAlves et al.ll2002t 
iMinniti et ail 1201 ll iGonzalez et alJl201 ll 120121) . Our compari- 
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Fig. 7. Left-hand panels: K s x (J - K s ) CMDs zoomed in the 
region around the Red Clump giants for areas #2, #5, #8 and #11, 
respectively, from the bottom to the top. The right-hand panels 
show the (J - H) x (H - K s ) CCDs for the same areas for all 
sources brighter than K s = 15 mag. 



son with models presented in Section 4 demonstrated that clump 
giants are the main population in the VVV bulge area, account- 
ing for up to 46% of the total sources at b — -3.0°. Even in the 
outer bulge, giants dominate at K s < 14 mag. 

Previous results have shown that the double RC in mag- 
nitude is the signature of the X-shaped MW bulge (e.g. , 
Mc Willia m & Zoccalil l20lot iNatafet all l20lot IZoccalil [2010; 

dto et al. 201 lbh . These analys es focused on inte rmediate 
Galactic latitudes (b < -6°), while ISaito etal.1 (1201 lbh mapped 
the X-shape structure for the entire bulge within 3.5° < \b\ < 
8.5°. The RC splits into two components for \b\ > 5, with the 
separation between the peaks increasing with the Galactic lati- 
tude, in agreement with the expectation of an X-shape structure 
seen almost edge-on, with the arms merging close to the Galactic 
centre. Our data are complementary to the previous results, and 
confirm the persistence of the X-shaped structure for b < -8°, 
as seen in the CMDs of the outer bulge (areas #1-3). At area 



#2 (b 9.0°) the double peak is seen very clearly in the VVV 

data, with the separation between peaks reaching AK S ~ 0.7 mag 
(see Fig. [7]i. 

Combining data from regions where the RC appears with two 
peaks at different separations produces an elongated structure in 
magnitude for the final dataset. On the other hand, the differ- 
ential reddening alo ng the bulge, with up to E(B - V) ~ 10 at 
the Galactic centre dGonzalez et al.l 120121) . causes the RC stars 
to spread along the reddening vector for more than A(J - K s ) > 
4 mag. The combination of these effects contributes to create the 
complex shape of the RC seen in our CMD for the whole bulge. 

5.2. The dwarf sequence 

The CMDs of the outermost bulge regions (areas #1-3) show at 
the right-most side a well-defined sequence of relatively faint, 
K s > 14 mag, and red objects (J - K s ) ~ 0.7 - 0.9 mag. The 
sequence merges with the disk sequence and the giant branch at 
lower magnitudes (K s > 16 mag), due to the large photometric 
errors in the bottom end of the CMD. The sequence is seen to be 
stronger towards negative longitudes, becoming more intense in 
the CMD for area #3. Interestingly, the sequence disappears at 
lower Galactic latitudes, and it is not visible in any other region 
of the VVV bulge. 

The comparison with models presented in Section 4 shows 
the red dwarf sequence as seen in the CMDs for the area at 
b — -9.5 deg. The models reveal that the sequence is composed 
of main-sequence stars belonging to the thin disk (see Fig. 0J. 
According to what is observed in CMDs for areas #4-6, the se- 
quence is not detected at b = -6.0°, because of the degeneracies 
introduced by crowding and differential reddening, which pro- 
duce the red dwarf sequence overlayed with the bulge RGB. 

We investigated the dwarf sequence in more detail and com - 
pared the colours of M-dwarfs provided bv lHewett et al.l J2006). 
Using colour transformations from UKIRT to the VISTA system 
and the mean extinction for the region, M-dwarfs from Ml to M6 
coincide in colour with the main sequence, while later objects 
are slightly redder. Fig. [4] shows these objects in the CMD, as- 
suming a distance modulus of (m^ s -M^ s ) = 7 mag (d ~ 250 pc). 

Interestingly, the red dwarf sequence is seen at K s > 14 mag, 
just beyond the limit of the 2MASS and DENIS photometry. 
This can explain why this was not noticed in pre vious near-IR 
studie s. Deep CMDs such as those presented by IZoccali et al.l 
d200l should reveal the sequence, but their analysis was fo- 
cused on fields at b ~ 6.0° where the red dwarf s equence is 
also a bsent in the VVV data, as described above. iLucas et ail 
(120081) reported a sequence with late K and M dwarfs in the 
UKIDSS GPS data, but at the opposite side of the Galaxy, at 
/ ~ 170 - 175 deg. However, the sequence is not clearly distin- 
guished from the RGB. The region at b =~ 9.0 - 10 deg seems 
to be ideal to survey these stars in the Galaxy, because the late 
dwarfs are seen in large numbers. The quality of the VVV pho- 
tometry allows one to easily isolate the sequence using colour 
criteria. K and M dwarfs are particularly important in the plan- 
etary transit searches during the variability campaign of VVV, 
because their small size allows a b etter contrast with transits of 
low-mass planets (ISaito et al.ll201 lal) . 



6. The CMD using reddening-free parameters 

The bimodal colour distribution of the red clump can be inter- 
preted as a reddening effect generated by the sampling of a large 
sky region with enough variation to produce the separation in 
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Fig. 8. KsXi J-Ks) CMD for area #09 (left-hand panel) compared 
with the CMD made using reddening-free indices (right panel). 
The notation is similar to that presented in Fig|3] 



two components, and it can be also created by the observation of 
two distinct populations seen at the same line of sight. 

In the VVV bulge area the RC spreading in colour, or even its 
splitting into a secondary peak, is seen across a large sky region, 
area #9 is the most drastic example. In this area the RC and even 
the RGB seems to be double in colour, with a separation of about 
A(J - K s ) ~ 0.2 mag. 

To test if the double RC is caused by differentia l reddening, 
we us ed reddening-free parameters provided by ICa telan et al. 
(1201 11) . based on the extinction law of Rieke & Lebofskv ( 1985) 
for the VISTA filters, namely m\ and C3, defined as 

m 4 = K s - 1.22(7-//), (1) 

c 3 = (J-H)-L47(H-K s ). (2) 

Figure [8]shows the K s x(J-K s ) CMD for area #9 (left panel) 
compared with CMD built using the reddening-free parameters 
TO4 x C3 for the same sources. It is very clear that the giant se- 
quence becomes single, and the RC is seen with a very sharp 
single peak. Not only the giant sequence changed, but all struc- 
tures are seen very concentrated (note the shorter colour range in 
the x-axis of the CMD built with the reddening-free parameters). 

7. Conclusions 

We presented the VVV colour-magnitude diagram of the 
Galactic bulge, the largest CMD ever published for a large ho- 
mogeneous data set, with 84M+ stars. The interpretation of the 
CMD yields important information about the MW bulge, show- 
ing the fingerprint of its structure and content. 

Stellar population synthesis models fit the data well, showing 
the predominance of main-sequence and giant stars in the outer 
bulge, which belong both to the bulge and halo. Thin- and thick- 
disk stars are also present in fewer numbers. In the inner bulge 
the CMD is dominated by bulge giants, which contribute up to 
46% of the total sources at b 3° 

The analysis of the outermost bulge area reveals a well- 
defined sequence of late K and M dwarfs, seen at (J - K s ) ~ 
0.7 - 0.9 mag and K s > 14 mag. These stars are particularly 
important in the planetary transit searches during the variability 
campaign of VVV. 

The RC appears to be double in magnitude in the outer bulge 
region due to the X-shaped structure of the MW bulge, with the 



separation between the peaks reaching AK S — 0.70 at b 9°. 

This result is complementary with previous analyses, and con- 
firms the persistence of the X-shaped structure for b < -8°. In 
contrast, in the inner part (b 3°) the RC appears to be spread- 
ing in colour, or even splitting into a secondary peak with a sep- 
aration of A(J - K s ) > 0.2 mag (area #9). The use of reddening- 
free parameters confirmed that this last feature is caused by red- 
dening effects. 

The CMDs of the Galactic centre are harder to interpret 
because they are strongly affected by reddening and extinc- 
tion. All populations are seen to be much redder, and the RC 
spreads along the direction of the reddening vector by A(J - 
K s ) > 2.0 mag. 
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